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ABSTRACT 

Space  craft  to  solar  system  bodies  must  deal  with 
multiple  deviations  under  the  in  uence  of  gravity.  The 
motion  of  satellite  in  one  of  the  types  of  conic  section 
strongly  depends  upon  escape  velocity.  The  goal  of 
this  paper  is  to  nd  the  trajectory  of  satellite  by 
analysing  orbit's  eccentricity  and  or-bital  parameters. 
This  paper  also  focuses  on  the  derivation  of  general 
orbit  equation  from  the  base  called  Newton's  law  of 
motion  and  Newton's  law  of  gravitation. 

Keywords:  Eccentricity,  Flight  path  angle,  True 
anomaly,  Perigee  and  Apogee  radius,  Trajectory. 

1.  INTRODUCTION 

Mechanics  is  a  branch  of  science  concerned  with  the 
relationship  between  force,  displacements,  energy  and 
their  e  ects  on  physical  bodies  on  their  environment. 
The  classical  mechanics  is  one  of  the  sub-  elds  of 
mechanics.  Classical  mechanics  is  concerned  with  the 
some  physical  laws  describing  the  movement  of 
bodies  under  the  in  uence  of  forces. 

2.  ORBITAL  MECHANICS 

Orbital  mechanics  is  also  known  as  astrodynamics.  It 
is  the  application  of  celestial  mechanics  concerning 
the  motion  of  rockets  and  other  spacecraft.  It  is 
concerned  with  motion  of  bodies  under  the  in  uence 
of  gravity,  including  both  spacecraft  and  natural 
astronomical  bodies  such  as  star,  planets,  moons  and 
comets. 
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2.1  TRAJECTORIES: 

Three  general  types  of  paths  are  possible  under  the 
gravitational  in  uence. 

An  elliptical  orbit  is  a  orbit  with  an  eccentricity  of  less 
than  1.  This  includes  the  special  case  called  circular 
orbit  and  its  eccentricity  is  equal  to  0. 

This  paper  examines  the  classifications,  cause;  types 
and  the  possible  solutions  of  conflict  in  organization. 
This  discovered  that  conflict  generates  considerable 
ambivalence  and  leave  many  practitioners  and 
scholars  quite  uncertain  about  how  best  to  cope  with 
it.  Conflicts  are  inevitable  in  human  life  and  also 
inevitable  in  organizations  even  between  countries. 
Conflict  occurs  in  organizations  as  a  result  of 
competition  for  supremacy,  scarcity  of  economic 
resource  and  leadership  style.  The  study  also  revealed 
that  conflict  in  organization  could  be  constructive  or 
destructive  which  can  lead  to  low  production  or  good 
solution  to  production. 

A  hyperbolic  trajectory  is  the  trajectory  of  any  body 
around  a  cen-tral  mass  with  more  enough  speed  to 
escape  the  central  object's  grav-itational  pull.  The 
orbital  eccentricity  will  be  greater  than  one 

A  parabolic  trajectory  is  an  orbit  whose  eccentricity 
will  be  equal  to  1.  A  body  travelling  along  an  escape 
orbit  will  coast  along  a  parabolic  trajectory  to  in  nity 
and  never  return. 
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2.1.1  ESCAPE  VELOCITY: 

Escape  velocity  is  the  velocity  required  at  a  given 
position  to  establish  a  parabolic  path.  Velocities 
greater  than  escape  velocity  result  in  hyperbolic 
orbits.  Lower  velocities  result  in  closed  elliptical 
orbits. 

2.2  ORBITAL  PARAMETERS: 

Apogee  (q)  refers  to  the  farthest  distance  between 
satellite  and  the  Earth. 

Perigee  (p)  refers  to  the  closest  distance  between 
satellite  and  the  Earth. 

3.  ORBIT  EQUATION 


Semi-Major  Axis  (a):  The  distance  from  the  center  of 
the  orbit  ellipse  to  satellite's  apogee  or  perigee  point. 

Semi-Minor  Axis  (b):  The  shortest  distance  from  the 
true  center  of  the  orbit  ellipse  to  the  orbit  path. 

Period  (T):  The  time  required  for  the  satellite  to  orbit 
the  Earth  once. 

Pitch  angle  or  ight  path  angle:  The  angle  between  the 
longitudi-nal  axis  (where  the  airplane  is  pointed)  and 
the  horizon. 


Two  objects  are  located  in  inertial  frame  of  reference  XYZ.  Each  of  the  bod-ies  is  acted  upon  by 
the  gravitational  attraction  of  the  other.  F 12  is  the  force  exerted  on  mi  by  m2,  and  F21  is  the  force 
exerted  on  m2  by  mi .  The  position  vector  Rq  of  the  center  of  mass  G  is  de  ned  by, 

_  mR+mR  ... 

Rg  —  11  2  2  !  (1) 

ml  +m2 


Absolute  velocity,  v 

G 


R0 


G 


ml  RIO  +m2R20  !  (2) 

ml+m2 


Absolute  acceleration,  a 


_  ^  _  miF^ao 
mi+m2 


!  (3) 


The  term  absolute  means  that  the  quantities  are  measured  relative  to  an  inertial  frame  of  reference. 
Let  r  be  position  vector  of  m2  relative  to  mi. 

r=R2  R!  !  (4) 

Let  ucr  be  the  unit  vector  pointing  from  mi  towards  m2, 

u  =  E  !  (5) 
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Cr  r 

where  r  =  krk  is  the  magnitude  of  r. 

The  force  exerted  on  m2  by  mi  is, 

F21  =  Gm1m^(  ub)  =  Gm^ub)  I  (fi) 

n  ■  wvv-  vwv-  ■  \  v  r 

j  j 

where  ucr  shows  that  the  force  vector  F2i  is  directed  from  m2  towards  mj. 

Applying  Newton's  second  law  of  motion  on  body  m2  is  F2i  =  m2R200,  where  R200  is  the  absolute 
acceleration  of  m2. 

=  m2Rm  |  (7) 

By  Newton's  third  law  of  motion,  Fi2  =  F2i,  so  for  mi  we  have, 

Gm,m.(uk)=m,R00  |(8) 

r2 

Equations  (7)and  (8)  are  the  equations  of  motion  of  two  bodies  in  inertial  space.  By 
adding  ,we  get  miRi°°  +  m2R2°°  =  0  !  (9). 

According  to  Equation(3) ,  the  acceleration  of  the  center  of  mass  G  of  the  system  of  two 
bodies  mi  and  m2  is  zero. 

Multiplying  Equation  (7)  by  mi  and  Equation  (8)  by  m2,  we  get 

Sn:n2  (jjLj  12  ,  ! 

12  =mnR  '  (10) 

Guii  m 

f  -  mi  m2  k  100  !  (11) 

Subtracting  the  equation(1 1)  from  (10)  yields 

m  m  (R00  R00}  =  Gm-mi  (rm_+_nr?luj_ !_( 1 2) 

I  I  2  ^ 

1  r* 

factor  mi  m2 

r"  =  G(m.+mJ  ur  !  (13) 

C 

Let  the  gravitational  parameter  be  de  ned  as 
=  G(m1  +  m2)  !  (14) 

Using  Equation  (14)  and  Equation  (5),  we  can  write  Equation  (13)  as 
r"  =  r3  r  !  (15) 
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This  is  the  second  order  di  erential  equation  that  governs  the  motion  of  m2  relative  to  m-|. 

Let  ri  and  r2  be  the  position  vectors  of  mi  and  m2  relative  to  the  center  of  mass  G.  The 
equation  of  motion  of  m2  relative  to  the  center  of  mass  is 

Gm,m;<V=nvr00  !(16) 

f  i 

where  r  is  the  position  vector  of  m2  relative  to  mi.  In  terms  of  ri  and  r2, 
r  =  r2  ri 

Since  the  position  vector  of  the  center  of  mass  relative  to  itself  is  zero,  we  get, 


in  1  r  1  +  m2r2  =  0 


ri 


—  m  2 


m  i  T2 


_  m  +m  m 
1  =  1  2  f2  1 


Substituting  this  back  into  Equation  (16)  and  using  (ucr)  =  r2=r2,  we  get 


Gmsim2  r=m  r°3 


[mi+m2)2r23  ^  ^ 


Simplifying, 


let  °=(  )3 

mi+m2 


m  1 


)3^r  =  r°°  (17) 


mi+m2  2  2  2 


r3°  r2  =  r002 


Similarly,  the  equation  of  motion  of  mi  relative  to  the  center  of  mass  is  found  to  be 


00  n  - 
-E3— .  ,r  "1  “ 


where  00  -  (  ) 


FTtl+fTli 


The  angular  momentum  of  body  m2  relative  to  mi  is  the  moment  of  ni2's  relative  linear  momentum  ni2r0 

1,  0 

H2=i  =  r  x  nri2r 


Divide  this  equation  by  rr^and  let  h  =  H2=i=m2,  so  that 

h  =  r  x  r°  !  (18) 

h  is  the  relative  angular  momentum  of  m2  per  unit  mass,  that  is,  the  sped  c  relative  angular 
momentum.  Taking  the  time  derivative  of  h  yields 
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=  r°  x  r°  +  r  f° 

But  r°  x  r°  =  0.  According  to  Equation  (15)  r"  =  r3  r,  we  have 

r  x  r00  =  r  x  ( r3  r)  =  r3  (r  x  r)  =  0 
=0  1(19} 
the  unit  vector, 

hb  =  — h  !  (20) 

vwv^  '  * 

Let  us  divide  the  relative  velocity  vector  r°  into  components  Vr  =  vrubr  and  Vp  =  vpubp  ,  where  ubr  and 
ubp  are  radial  and  perpendicular  unit  vector 

b 

h  =  rubr  x  (vrubr  +  VpUbp)  =  rvph  h  =  rvp  !  (21 ) 

The  angular  momentum  depends  only  on  the  perpendicular  component  of  the  relative  velocity. 

During  the  di  erential  time  interval  dt  the  position  vector  r  sweeps  out  an  area  dA. 

From  the  gure  the  triangular  area  dA  is  given  by 


1  1  1 
dA  =  _2  x  base  x  height  =  _2  x  vdt  x  r  sin  =  ~2  r(vsin  dt)=1 2rvpdt 

Therefore,  using  Equation  (21)  we  get 

dA  _  h 

dt  -  2  ■  (22) 

The  vector  identity  also  known  as  the  bac  -  cab  rule: 

A  A  x  (B  x  C)=  B(A:C)  C(A:B)  !  (23)  r:r  =  ?  !  (24) 

dt  (r:r)  =  2r-dt 
But 

d  ,  v  _  d  r  d  r  dr 

dt  (r.r)  -  r  dt  +  dt  r  -  2r^ dt 
Thus,  we  obtain  the  identity 
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r:r  =  r:r  !  (25a) 

Since  r°  =  v  and  r  =  krk  ,  this  can  be  written  alternatively  as 

kdr  k 

r:v  =  krk  dt  !  (25b) 

Let  us  take  the  cross  product  on  both  sides  of  Equation  (1 5)  r"  =  r3  r  with  the  speci  c 
angular  momentum  h:  — 


r"  x  h=  r3  r  x  h 


!  (26) 


Since  dtd  (r°  x  h)  =  r°°  x  h  +  r°  x  h°,  the  left-hand  side  can  be  written 


°o  ,  d  .  o  ,  .  o  .  o  . 
r  x  h  =  dt  (r  x  h)  -  r  x  h  !  (27) 


But  according  to  Equation  (19),  the  angular  momentum  is  constant  (  h  =0),  so 

00  d  o  M 
r  x  n  —  dt  (r  x  nj 

The  right-hand  side  of  Equation  (26)  can  be  transformed  into: 

r13_rxh  =  r13_(rx(rxr°)) 


But 


Therefore 


=  r  j  [r(r:r°) 

r°(r:r)] 

=  f1 3  [r(rr°) 

0/  2vi 

r  (r  )] 

(rr° 

r°  r) 

(r°r  rr°) 

rr° 

)  =  r2 

=  ( 

,  d 

i  r  x  h  =  dt 

(”r  )  !  ( 

r2 


Substituting  Equations  (27)  and  (28)  into  Equation  (26),  we  get 

d£(r°*h>  =  dtS{  Sj 
^Ljr0  x  h  -  :)  =0 


Taking  the  dot  product  of  h  on  both  sides  of  Equation  (29)  yields 
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(r°xh)dL 


r.h 

(  V  )  =  C  .  h 


Since  r'  x  h  is  perpendicular  to  both  r'  and  h,  so  that  (r'  x  h)  .  h=  0.  Likewise,  since  h  =  r  x 
r'  is  perpendicular  to  both  r  and  r',  it  follows  that  r .  h  =  0.  Therefore,  we  have  C  .  h=  0, 
i.e.,  C  is  perpendicular  to  h,  which  is  normal  to  the  orbital  plane. 

Rearranging  Equation  (29) 


r 

r  +  e  = 


1  (r°x  h)  !  (30) 


where  e  =  C=  .  The  vector  e  is  called  the  eccentricity.  The  line  de  ned  by  the  e  is 
commonly  called  the  apse  line,  o  obtain  a  scalar  equation,  let  us  take  the  dot  product  of 
both  sides  of  Equation  (30)  with  r 

r  r 

_  r:  +  e  =  1  (r:(r°x  h  ))  1(31) 


we  know  the  identity 


to  obtain 


A:(BxC)  =  (AxB):C 


r:(r°xh)  =  (rxr°):h  =  h:h  =  h2  !  (32) 

Substituting  this  expression  into  the  right-hand  side  of  Equation  (31),  and  substituting  r . 
r=  r  on  the  left  yields 

r+r:e  =  h2  1(33) 


From  the  de  nition  of  dot  product  we  have 

r:e  =  re  cos 

In  terms  of  the  eccentricity  and  the  true  anomaly,  we  may  write  Equation  (33)  as 

tl2 

r  +  re  cos  = 


1  +e  cos 

This  is  the  orbit  equation,  and  it  de  nes  the  path  of  the  body  m2  around  mi,  relative  to  mi. 

3.1  Flight  path  angle: 

The  ight  path  angle  is  the  angle  that  the  velocity  vector  v  =  r'  makes  with  the  normal  to 
the  position  vector. 

.  V 

tan  = _ [ 

Vp 
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tan  = 


1  +e  cos 


3.2  Period  of  orbit  (ellipse): 


At  periapsis,  =  0 


&=  T.7-  '  P4) 


At  apoapsis,  =180 


la  "  !  (35) 

WvW 


Let  2a  be  the  distance  between  from  periapsis  P  to  apoapsis  A,  2a  = 

A. 

2a  =  —  —  4*-^ 

1+e  1  e 


a  = 


1 


1  e 


Area  of  ellipse  =  A  =  ab 
By  Kepler's  second  law, 


b  =  i  e2 


dA  h 
dt“  2 


A  -  2  t 


where 


A  =  ab  and  t  =  T 


Substituting  the  values  of  a  and  b, 


ab  -  2  (T  ) 


T  =2ab 


T  —  2 —  /P .  = 

I  -  2  (  1  e2  ) 


2  3=2 


T  =  pa 


From  the  equation  of  semi  major  axis,  we  have  h  =  a(1  e2) 


IJTSRD  |  Available  Online  @  www.ijtsrd.com  |  Volume  -  2  |  Issue  - : 


rp  +  ra 


Jan-Feb  2018 


Page:  1270 


International  Journal  of  Trend  in  Scientific  Research  and  Development  (IJTSRD)  ISSN:  2456-6470 


3.3  Eccentricity: 

From  equations  (34)  and  (35)  we  have, 


rp  —  1  e 

WvVVv'  —  — 

<A  1+e 

0  t‘.  V?V. 


3.4  NUMERICAL  ILLUSTRATIONS: 

1 .  A  satellite  orbits  the  earth  with  a  perigee  radius  of  7000  km  and  apogee  radius  of  70000  km. 

Calculate  eccentricity  and  period  of  orbit. 

Solution: 

e  =  Up 

_ Yv  Yv 

ISrfSB- 

=  7DDDD  7DDD 
70000+7000 

e  =  0:B1B2 

Since  e  1,  it  will  be  an  ellipse. 


we  know  that  r  = 

1  +e  cos 


1+e 


7000  =  — - 

398600  1+0:8182 

h  =  71226  km2/s. 

Period  of  orbit,  T  =  2  2  (Pjh  e_2  )3 


T  —  2  71226  ^3 

3986002  1  0:81 822 

T  =  751 1  seconds 
T  =  20:89  hours 

2.  At  two  points  on  a  geocentric  orbit  the  altitude  and  anomaly  are  Zi  =  1700  km  ,  1  = 
130°  and  Z2  =  500  km  ,  2  =  50°.  Find 

a)  eccentricity 

b)  angular  momentum 

c)  perigee  radius 
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d)  perigee  altitude 

e)  semi  major  axis 

Solution: 

ri  =  Re  +  zi  =  6378  +  1700  =  8078  !  (1) 
r2  =  Re  +  z2  =  6378  +  500  =  6878  !  (2) 

h2  i 

we  know  that  r  = - 

1  +e  cos 

h2  1 

(1)  )  8078  =  (398600)  1+ecos  130° 

h2  =  3:219  xIO9  2:069  x  109e  !  (3) 

h2  1 

(2)  )  6878  =  (398600)  1+ecos  50° 

h2  =  2:742  x  109  +  1:762  x  109e  !  (4) 

by  (3)  and  (4),  3:219  x  109  2:069  x  109e  =  2:742  x  109  +  1:762  x  109e 

3:219  2:742  =  e(1:762x  2:069) 
e  =  0:477=3:831 
e  =  0:1245 

Since  e  1,  it  will  be  an  ellipse. 
h2  =  3:219  x  109  2:069  x  109(0:1245) 
h2  =  2961409500 
h  =  54418  km2 /s 


2 

54418  1 

rP  “  398600  1+0:1245 
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rp  =  6607  km 

d)  zp  =  rp  Re 

zp  =  6607  6378 

zp  =  229  km 
h2  1 

e)  ra  =  — 

1  e 


_  _  544182  1 

<  a  - - 

398600  1  0:1245 


ra  =  8486 

g  =  ra+rp  2 


g  =  8486+6607  2 

a  =  7546:5 

3.  At  a  given  point  of  spacecraft's  geocentric  trajectory,  the  perigee  alti-tude  is  300  km, 
the  speed  is  15  km/s  and  the  ight  path  angle  is  50°  Show  that  path  is  hyperbola  and 
calculate 

a)  hyperbolic  excess  speed 

b)  angular  momentum 

c)  true  anomaly 

d)  eccentricity 

Solution: 

a)  rp  =  Re  +  zp  =  6378  +  300  —  6678 


vesc  -  r 

2(398600) 

_q - 

veSC  “  6678 

vesc  =  10:925 

2  2 
Vinf2  -v  -  ves  c 

=  225  119:377 


@  IJTSRD  |  AvailableOnline@www.ijtsrd.coml  Volume -2  |  Issue -2  |Jan-Feb2018 


Page:  1273 


International  Journal  of  Trend  in  Scientific  Research  and  Development  (IJTSRD)  ISSN:  2456-6470 

inf2 =105:623 

Vinf=1 0:277  km2=S2 

b) vr  =  rsin  =  15  sin  50°  =  11:49 

vp  =  rcos  =  15  cos  50°  =  9:641 

h  =  rvp  =  6678(9:461) 
h  =  64383  km2/s 

c)  we  know  that  r  =  —  1 

1  +e  cos 

6678  =  643832  1 

398600  1+e  cos 

e  cos  =0:5573  !  (1) 
v  =  e  sin 

r  h 

398600 

1 1 :49  =  64393  e  sin 

e  sin  =  1 :856  !  (2) 

,  e  sin 

tan  = - 

e  cos 

tan  =„v5573856 

=  73:28° 

d)  (1)  )e  cos  73:28°  =  0:5573 
e  =  1:9379 

Since  e  1 ,  it  will  be  hyperbolic  orbit. 

(5).  The  perigee  of  a  satellite  in  a  parabolic  trajectory  is  9000  km.  Find  the  distance  d 
between  points  Pi  and  P2  on  orbit  which  are  11000km  and  22000  km  from  the  center  of 
earth. 

ho  1 

Solution:  We  know  that,  rP  = 

1+e 
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Figure  1:  parabola 

Since  it  is  a  parabolic  trajectory  its  eccentricity  will  be  one. 

rp,=^  _ 

1+1 

h  =  p2V 

£  - 

h  =  2(398600^9000) 

h  =  84704km2/s 

True  anomalies  of  points  Pi  and  P2  are  found  using  orbit  equation, 

we  know  that  r= 

iiooo  =  §™1 _ 1_ 

398600 1  +cos  -i 

1  =  50:48° 

22000  =  _ 1_ 

398600 1 +C05  2 

2=  100:48° 

-7  1  =  100:48°  50:48° 

VWVTO  u 

=  50° 

From  law  of  cosines  from  trignometry, 
d2  =  1 1 0002  +  220002  2(1 1 000)(22000)  cos 

50°  d  =  17143  km 
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